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ABSTRACT 


The present study takes the initial steps of establishing a better flame trench 
design criteria for future National Launch System vehicles. A three-dimensional 
finite element computer model for predicting the transient thermal and structural 
behavior of the flame trench walls was developed using both I-DEAS and 
MSC/NASTRAN software packages. The results of JANNAF Standardized Plume Flowfield 
calculations of sea- level exhaust plumes of SSME, STME and ASRM were analyzed for 
different axial distances. The results of sample calculations, using the 
developed finite element model, are included. The further suggestions are also 
reported for enhancing the overall analysis of the flame trench model. 
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SUMMARY 

Pads 39A and B ate being studied “ 

can be launched without damaging t e exi family of vehicles which are 

designated as NLS (National Launrf . System^ NLS^is^a family^ ^ igurations . The 

built around a common core ^ ^ determined ye t. Options include 

engines for the core an csSME) modified Saturn V engines (F-1A) , Russian 

using old shuttle main engi C ’ fla me trench capacity will be an 

t^Sed; spacebar. Due to the importance of this problem 
an accurate, reliable solution is required. 

The purpose of this summer research stud y ^ 

the existing flame trench and to esta is problem requires the 

for future NLS launch vehicles. The solution or y 

following steps: 

1 A three-dimensional finite-element model of the flame trench walls, 
o' The exit flow conditions of the rocket engine 

Z. me exu j-j- . f rha 5P p.1evel plume characteristics. 

3. The complete analysis of the se P d pressure rates on the 

4. The plume impingement code to predict the heating an p 

flame trench. 

The above steps are carried out as far as possibU in ^ ^i^time 

available. As a first step three '< and structural behavior of the 

developed for pred ^ C p ^| and MSC/NASTRAN software packages . Secondly, the JANNAF 
flame trench using I- DEAS and MSC/ ^ obtain the gas dynamic structure 

Standardized Plume f : lowf] ield < ( SPF) rocket exhaust flow fields. The analysis 

of sea- level exhaust plumes of differ heating rates and pressure rates 

is also needed to gain a better un ers described above are not completed 

on the flame trench walls. Since th P this report should be looked 

fully at the present time, the results’ Results will be used later to improve 
at in a qualitative manner. However, these resuits will for improving 

the required solution. ^ ««rth« S ^ concluded that the three- 
the overall analysis o e the fi ame trench can be tested for any given 

modifications explained in section V. 
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I. INTRODUCTION 

The existing flame trench at Launch Pad 39A was built in 1967 f ° r 
space program. An identical launch pad was bui Itlnl * trench 
pads were converted to be used with the space «huttl« tinWl.JL* ^ 

has not been modified since ^ its ^ be launched wit hout damaging the 

to determine if a new fami y NLS (National Launch 

existing structure. The new vehic are built around a common core using 

System) . NLS is a family of vehicles which are ^ ^ and boosters 

different external booster Jjonfigu rat • f d shuttle main engines (SSME) , 

have not been determined. Options includ^^g^ ^ > and solid fuel 

roosters. The" flame S capacity’ will be an important factor in the decision 
process . 

The flame trench is built •* “1 S^etdeep*. rtlTTsl 

vehicle te divert the flow t. the 

north and south exists from the flame trenc . 

The flame trench is limited in ^ 

« of crack and excessive amounts 

of temperatures could damage the refractory coating. 

KSC-STD-Z-0012B , Standard for Flame Deflectc orjtesi 
criteria standard currently aval a e o ^ heig ht, width and distance away 
establishes minimum design re qui remen comply with. Currently, the 

£ 

trench must then be sized to accommodate the flame 

deflector . 

The purpose o£ this study is to investigate the 
flame trench for different exhaust t°nft gut at ions an sCu dy is carried out 

trench design criteria ^ 1 shows the general overview of 

based on the drawing number ana rxg 

the PAD39B. 
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Main thrust of the summer research study can be summarized as: 

(1) Survey literature of available computer codes for the problem solution 
deluding y the noz,U design, plume flow, plume impingement and finite-element 

programs . 

(2) Developing a three-dimensional finite-element model using I-DEAS and 
MSC/NASTRAN softwares for the flame trench walls. 

(3) Using Joint Army. Navy, NASA, Ait Force (J MUF V 
code (SPF) for predicting the gas dynamic structur 

plumes . 

(4, Better prediction o£ the pressure and heating rates on the flame trench walls 
using Source Flow Plume Impingement Program (SFPLIMP). 
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II. descriptions of the computer codes required for the study 

The computer codes required to solve the problem are briefly described 

below: 

2.1 RAMP CODE 

The Reacting multiphase (RAMP) code iSTvMUby ”of ^roMemt'assocIaild 
under government funding and use to so > ve code is capable of 

„ith real gas, supersonic, compressible flou (Ref. 1 . nterface „ ith che 

performing rocket nozzle flow ca cu a 10 ^ considers the exchange of 

JANNAF CODE. The code is fully coup ' ® t ™ie plumes. The results of the nozzle 
momentum and energy between the gas P radial distributions of 

proper tie^at^the e^ pTane S i^it^l Editions for the JANNAF code. 

2.2 JANNAF CODE 

There has been remarkable P*°**®®®.*” ^ V ®^ P ^^ at p henomena in the past 
of the complicated structure an ut;er programs are available to give the 

decade. A number of review papers state-of-the-art in rocket exhaust plume 

researcher a general overview of the state 

analysis (Refs. 2 through 4). 

r- ( t akima'F'i Standardized Plume Flowfield 

Joint Army, Navy, NASA Air °*ce (JA^ and industr y standard 

(SPF) is a computer code which is a p ° hase low altitude rocket exhaust 

for predicting the gas dynamic s true tu P flowfield regions. The 

plumes. Figure 2 illustrates a sch ?““%°/ te “ e v J rsion (SPF-II) of the JANNAF 
computational methodology employed m the latest v 

code is described in Ref. 5. 

The code main!, contains 

Capturing and Turbulent Mixing compo ^ initial conditions and external 

input data (run P arameterS ’ (JANNAF thermodynamic and chemical kinetic data 

data a “ bedded " * aubroucine ' 

savor n f the code provides an inviscid, frozen 
The Shock-Capturing Component of t P It uses shock -capturing to 

“ r Lr S " 

approximately £ three times the distance to 

the Mach disc. 

v p CPF nrovides a viscous, turbulent mixing 
The Turbulent Mixing Component of SPI ? chemistry. Gas/particle 

solution of the plume shear layer f viscosity and thermal 

m ° delinS thr ° USh "" 

component . 
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2.3 SFPLIMP CODE 

The Source Flow Plume ImpingementJrogrM 
to compute the forces , moments an i e g ig capable of modeling on a wide 

bodies Immersed rnto P lu “ fl ° lume a „d processing dace, user must 

rectangles a'nd^.T-d 'ensTonal shapes (circular cones, circular 

cylinders, spheres and polynomials of revolution). 

2.4 I-DEAS FEM PACKAGE 

I-DEAS Finite Element Modeline (FEM) a \ l ^ S t ^ tl “ a er 1 '° ds bUl a 1 nd boundary 

element model including ^“/^“^"^pTe-processing ' and post-processing 
conditions. This package _ DEA c Model Solution. The model 

features. Solving the model .generally runs in ^ D “ S ^1 . trans£er . 

solution can solve problems in linear ^ S J an sfer. I-DEAS Model 

SUl^tsld, solver package. 

2.5 MSC/NASTRAN CODE 

MSC/NASTRAN is a large scale g enera lP^P°s e ^^ i^ t e°element P method. This 
solves a wide variety of engineering P*° ie * t heat transfer with constant or 

code can analyze linear, non-linear vg boundary conditions. The code can 

temperature -dependent convective dependent thermal conductivity material 

^.r^irr 1 ^: — - — - 

in transient solution algorithm. 

MSC/MASTEAM operates in a hatch mode . Input 

in the form of a card deck or “^ntrol dec^ tS^ purpose of this deck is to 
preparing the cards for •*= »!•<> the 

identify the job and the type of . , exeC uted, such as maximum time 

general conditions under which the job is to b “ecu': ^ ^ the data 

allowed, type of system defines the load case, selects the 

cards for the case control dBC^ Th ^ requests for plotting. The Bulk 

Deck % of organizing £2 "yt 

ass ££.2r££ s £ 
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III. METHODOLOGY 

Figure 3 shows the computational scheme required for the problem solution. 
The JANNAF computer code requires the exit conditions of the engine nozzle to 
start the calculations for sea-level plume definition of the proposed exhau 
configuration. The results from the JANNAF code must be run to Predict; th 
loading conditions of the finite element model. The I-DEAS and MSC/NAST co es 
must be then used to determine the capacity of the flame trench for different NLS 

configurations . 

3.1 FINITE ELEMENT MODEL 

The flame trench facing to the north side of Pad 39B is modeled by a three- 
dimensional finite element model starting at the base of the exhaust deflect • 
running approximately 50 feet toward the north along the centerline of the plum 
footprint Figure 4 shows the three-dimensional model used in the study. The 
model is built of from 280 nodes and 162 solid brick type of elements^ The model 
uses physical and thermal material properties of concrete and brick The concre « 
stiffeners shown in drawing 79K10338 were modeled with the restraint set as shown 

in Figure 5 . 

3.2. HEATING LOADS AND PRESSURE LOADS 

The impingement of the Space Shuttle Main Engines (SSME) and Solid Rocket 
Motors^fSRJL) or Advanced Solid Rocket Motors (ASRM) gaseous exhaust plumes on 
Launch Complex 39 (LC-39) results in pressure loads and convective heating 

Serous components. Reference 6 presents the results of a study to determine the 
thermal and pressure environment imposed on elements of LC-39 by impingemen 
the Space Shuttle Vehicle (SSV) exhaust plumes at the sea level for a maximum 
north drift trajectory. The plume properties given in this reference are us 
establish pressures and heating loads to objects immersed in the ASRM or SSMt 
exhaust plume. Reference 6 do not take into account for radiation from pi 
exhaust turned by the Mobile Launcher Platform (MLP) and recirculation. 
Prediction of this water-cooled launch pad radiation is not readily computed due 
to tie extreme changes of rocket plume flow field properties and recirculated 
flow field environments during launch. 

Both SRM and SSME flame trenches are subjected to heating due to vehicle 
exhaust gases that have flowed down the trench after having impinged on the flame 
or Arigorous calculation of the thermal and pressure environments of 
obfects subjected to indirect impingement requires a three-dimensional analysis 
of the flow J field which would be far beyond this study due to time limitation^ 
Therefore° a simplified approach is used for the heating and pressure loads 
required for the study. The heating rates were calculated for each point on 
flame trench as described in Reference 6, assuming that the heating rate at th 
point where the plume initially strikes the flame deflector varies with the 
distance down the flame deflector and flame trench, taken to the 0.2 power. 
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The heat loading conditions, used in the present study, lnduda the face 
heat influx distribution to fifty-four of the elements as Presented in T 1^1 
as explained in the previous paragraph and face convection applie 
"rfaS of th" model. The heat transfer coefficient for face convection was 
calculated based on the recommendation of Churchill and Chu ( . ) • 

The results of the gaseous and particle plume impingement pressure loads 
from the piles are also included in Reference 6. It is assumed that the gaseous 
flow and particles travel parallel to the trench walls. The estimates of pressure 

loads to the flame trench are presented as a result of ind ^ r ® C ^ ^"Jeffrenc e 6 
face pressure level of 2 atmospheres is used as recommended in Reference 6. 
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IV. RESULTS 


presented in this section should be looked at qualita iv y. 

Figure 6 depicts the 

pressure loading of 2 atmospheres for the static run y splacelI ent 

indicates the maximum displacement of 002 inches. stiffeners at the 

contour levels ^Following* program execution the deformed shape of the 

ifshrln Figure 7* ft U STIt^ 

manner which agrees with Its expected deformation, tlgur 
contour levels for the same loading condition. 

Pigure 9 illustrates the steady-state resulting 

model based on the applied face he ..hair were observed at alumina brick 

Maximum temperatures MSC/NASTRAM 

close to the deflector p ^ : Table 2 presents the results of 

t^Vransielt 1 finite element program during the first 

three seconds of launch. 

Both I-DEAS and MSC/NASTRAN are capable of computing thermal^ stress analyst 

due to the thermal loading conditions. n ,. .. uti on i n the model due to the 

analysis) will calculate the temperature di^rilmti™^ "„ ctural analysis) 
given thermal loads and boundary co . Ioada co a structural analysis 

uses these calculated temperatures as temper.ture^loads Flgur e 

S an^l^illustrate^the^displaceraent and stress contour lines, respectively. 

The sea level plume definitions were c ^“ la “ d Xsbu.Tl ttaotghT'inXidf 
an exit plane start line generated with the RAM pressure) of ASRM, SSME and 

feet . 
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V. CONCLUSIONS AND RECOMMENDATIONS 

* i zs e= d c^ g p r^ s y r^ p — ofj 

flame trench capacity of diff h use a effectively to analyze 

finite element model presented in this s u y trenc h for proposed loading 

conditions^ Further^nprov^ents on^the FE^^odel may include: a<Mition of ^the 

the resnlts with a 

better mesh refinement of the FEA model. 

Suggestions for completion of the problem are listed below; 

- Use of reacting multiphase nozzle code for Saturn V and Russian RD-170 engines. 

. Run of the JANNAF Plume code for Saturn V and Russian RD 170 engine exit nozzle 
conditions . 

' Sdift^Tea^U'lo"^ 

- Setting temperature and stress limits for the flame trench model. 
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appendix 


ID DIRECT TRANSIENT, P3DM0D1 

$ 

$ MAXIMUM CPU TIME ALLOWED FOR THE JOB 

$ 

TIME 10 $ UNITS ARE IN MINUTES 

I THE THERMAL ANALYZER PORTION OF MSC/NASTRAN IS TO BE USED 

$ 

APP HEAT 

$ 

$ THE TRANSIENT ALGORITHM IS TO BE USED 

$ 

SOL 89 

CEND $ END OF EXECUTIVE CONTROL DECK 

$ 

$ END OF EXECUTIVE CONTROL . . . START CASE CONTROL 
$*■*"**** 

$^ TL REQUEST I SORTED°aMD^ UNSORTED CArITeS'o-S to SEE COMMENTS 
! r™iS CARD IS OMITTED ONLY THE SORTED BULK DATA WILL APPEAR 
$ COMMENT CARDS WOULD THEN NOT BE PRINTED 
ECHO-BOTH 

! SPECIFY THE SUPERELEMENTS TO BE RUN. THIS CARD IS REQUIRED 

$ $ AND^ SPECIFIES ALL FOR THIS MODEL THERE IS BUT A SINGLE 

$ SUPERELEMENT WHICH ENCOMPASSES THE ENTIRE MODEL 

SEALL=ALL 

DLOAD=300 

TEMP (ESTI) =400 

IC=400 

TSTEP=500 

$ SELECT THE OUTPUT DESIRED 

$ OUTPUT 

THERMAL-ALL 

SUPER-ALL 

FLUX-ALL 

SUBTITLE-CASE SET 2, LOAD SET 2 
$ END CASE CONTROL DECK 

***** BULK DATA RECEIVED FROM THE I-DEAS ****** 

***** 1085 LINES OF ENTRY ****** 

TLOAD2 , 300 , 2 , , 0 . ,l.+6,0. ,0. ,+TLl 
+TL1 , 0 . , 0 . 

TEMPD.400,70. 

TSTEP , 500 , 10 , .2,1 
PARAM AUTOS PC YES 
PARAMPOST -2 
ENDDATA 
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TABLE 1 

The Face Heat Influx Distribution on the 
Brick Elements using I-DEAS software 


HEATING RATE (Ref. 6) 
BTU/(FEET**2 SEC) 


DISTANCE FROM THE 
DEFLECTOR BASE 
(FEET) 



3.5 

10.5 

2 

240.16 

282.44 

6 

192.81 

226.76 

10 

174.04 

204.69 

15 

160.55 

188.83 

21 

150.12 

176.57 

27 

142.70 

167.83 

33 

137.08 

161.22 

40.5 

131.91 

155.14 

47.5 

128.86 

151.54 


DISTANCE FROM 
THE FLOOR 
(FEET) 


12 

.5 

24 

.5 

31. 

5 

38.5 


372. 

68 

372. 

68 

282. 

44 

240. 

56 

299. 

20 

299. 

20 

226. 

76 

192. 

81 

270. 

09 

270. 

09 

204. 

69 

174. 

04 

249. 

16 

249. 

16 

188. 

.83 

160. 

55 

232. 

97 

232. 

97 

176 

.57 

150. 

12 

221. 

45 

221. 

45 

167 

.83 

142. 

.70 

212. 

73 

212. 

73 

161 

.22 

137 

.08 

204. 

.71 

204. 

.71 

155 

.14 

131 

.91 

199. 

.96 

199. 

.96 

151 

.54 

128 

.86 
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TABLE 2 


Results of the Transient Analysis with 
MSC/NASTRAN Program 

TEMPERATURE (DEGREES IN FAHRENHEIT) 
DISTANCE FROM THE DEFLECTOR BASE 


TIME 

(SEC) 

X=0' 

X=12 ' 

X=24 ' 

X=36 ' 

X*=50 ' 

0.0 

100 

100 

100 

100 

100 

0.2 

227 

185 

174 

168 

165 

0.4 

505 

371 

337 

318 

309 

0.6 

731 

523 

471 

441 

426 

0.8 

893 

634 

568 

530 

511 

1.0 

1037 

734 

655 

610 

588 

1.2 

1160 

819 

730 

679 

654 

1.4 

1272 

898 

800 

743 

713 

1.6 

1376 

972 

864 

802 

772 

1.8 

1474 

1043 

926 

859 

826 

2.0 

1568 

1111 

985 

914 

878 

2.2 

1659 

1178 

1043 

967 

929 

2.4 

1747 

1243 

1100 

1019 

979 

2.6 

1834 

1308 

1157 

1071 

1029 

2.8 

1920 

1373 

1212 

1122 

1078 

3.0 

2003 

1436 

1268 

1173 

1126 
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TABLE 3 

Sea Level ASRM , SSME and STME Exhaust Plumes 
Axial Mach Number Distributions 


Axial Distance 
(FEET) 

50 

100 

150 

200 

250 

300 

350 

400 

450 

500 

550 

600 


MACH NUMBER 


ASRM 

SSME 


2. 

363 

1 . 

466 

2. 

363 

1 . 

145 

2. 

358 

0 . 

615 

2. 

096 

0 . 

387 

1 . 

.758 

0 . 

.278 

1 . 

.467 

0 . 

.216 

1 

.224 

0 

.177 

1 

.026 

0 

.150 

0 

.869 

0 

.131 

0 

.743 

0 

.116 

0 

.642 

0 

.105 

0 

.561 

0 

.096 


STME 


2.690 

1.441 

0.812 

0.494 

0.344 

0.262 

0.211 

0.177 

0.153 

0.134 

0.120 

0.109 


TABLE 4 

Sea Level ASRM, SSME and STME Exhaust Plumes 
Axial Temperature Distributions 


Axial Distance 
(FEET) 

50 

100 

150 

200 

250 

300 

350 

400 

450 

500 

550 

600 


TEMPERATURE 
(DEGREES RANKINE) 

ASRM 


3879 

3879 

3898 

4233 

4607 

4848 

4973 

4988 

4904 

4740 

4496 

4175 


SSME 

STME 

4918 

4060 

4564 

4849 

3198 

4194 

2134 

2890 

1597 

2080 

1308 

1635 

1134 

1371 

1021 

1202 

943 

1086 

885 

1003 

842 

941 

808 

893 
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TABLE 5 


Sea Level ASRM, SSME and STME Exhaust Plumes 
Axial Pitot Pressure Distributions 


PITOT PRESSURE 
(PSIA) 


Axial Distance 
(FEET) 

50 

100 

150 

200 

250 

300 

350 

400 

450 

500 

550 

600 


ASRM 

SSME 

105.8 

51.9 

105.8 

31.4 

105.3 

18.6 

84.5 

16.2 

61.6 

15.5 

45.4 

15.2 

34.3 

15.0 

27.3 

14.9 

23.1 

14.9 

20.5 

14.8 

18.9 

14.8 

17.9 

14.8 


STME 


97.6 

43.2 
21.8 

17.2 

15.9 
15.4 

15.2 
15.0 

14.9 
14.9 
14.9 
14.8 
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gure 2 A Schematic Representation of the Plume Flowfield Regions 
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Figure 3 


Computational 


Scheme 
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Figure 


6 


Displacement Contour Lines 


of Constant Pressure Loading 
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Figure 9 Steady-State Temperature Contours of the FEA Model 
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